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ABSTRACT 

We analyze the interaction of a radiation-dominated jet and its surroundings using the equations 
of radiation hydrodynamics in the viscous limit. In a previous paper we considered the two-stream 
scenario, which treats the jet and its surroundings as distinct media interacting through radiation 
viscous forces. Here we present an alternative boundary layer model, known as the free-streaming jet 
model - where a narrow stream of fluid is injected into a static medium - and present solutions where 
the flow is ultrarelativistic and the boundary layer is dominated by radiation. It is shown that these 
jets entrain material from their surroundings and that their cores have a lower density of scatterers 
and a harder spectrum of photons, leading to observational consequences for lines of sight that look 
“down the barrel of the jet.” These jetted outflow models may be applicable to the jets produced 
during long gamma-ray bursts and super-Eddington phases of tidal disruption events. 

Subject headings: galaxies: jets - gamma-ray bursts: general - radiation: dynamics -- relativistic 
processes 


1. INTRODUCTION 

Particles, magnetic fields, and radiation all contribute 
to the propulsion of relativistic jets, though the relative 
contribution of each is still an open matter of debate. In 
certain situations, however, the mechanism responsible 
for launching the jet operates simultaneously with the 
release of a large amount of energy in the form of radi¬ 
ation, making it likely that photons dominate the bulk 
energetics. This scenario occurs, for example, during the 
super-Eddington phase of jetted ti dal disruption events 
(TDEs), such a s Swift J1644+57 (|Burrow/eLalJ _l20j4l 
Zauderer et al.l[2Mlh and Swift J2058+05 ( Cenko et al.l 


2012f) . Radiation-dominated jets should also be present 


in the collapsar picture of long gamma-ray bursts (GRBs; 
Rees fc Meszaroslll992t iWooslevlIlQQ.Jl iMeszaros fc ReesI 


1993t iMacFadven fc WooslevI 119991 iPiranI 1200411 . where 
the energy released in the form of gamma-rays is ul¬ 
timately derived from accretion onto a black hole, the 
associated accretion luminosity exceeding the Edding¬ 
ton limit of the hole by more than ten orders of mag¬ 
nitude. In both of these cases, the propagation of 
the radiation-dominated jet is modulated by the pres¬ 
ence of a radiation pressure-supported environment; 
for super-Eddington TDEs, this environment is in the 
form of a hi ghly inflated, quasi-spherical torus of fall¬ 
back debris (jCoughlin fc Begelmanll2014afl. while a “co¬ 
coon” of shocked jet rn aterial (|Morsonv et al.l l2007t 
iLopez-Camara e t al.l201 3D and the overlying stella r enve- 


lope itself (|Matzned[2003HWooslev fc Hegenl2006ll serve 

as the confining medium fo r GRBs. _ 

In a companion paper (|Coughlin fc BegelmanI 120151 
hereafter Paper I), we presented a model that describes 
the viscous interaction of a radiation-dominated, rela¬ 
tivistic jet with its surrounding medium. In that analysis 
we treated the jet and its surroundings as two separate 
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fluids, interacting with one another via small anisotropies 
in the comoving radiation field that are explicitly ac¬ 
counted for in the equations of radiation hydrodynam¬ 
ics in the viscous limit. This two-stream approxima- 
tion, in agreement with the non-relativistic analysis of 
lArav fc BegelmanI (|1992f ). demonstrates the manner in 
which the shear between the two fluids carves out a re¬ 
gion of low density material within the boundary layer 
between them. We also deduced the dependence of the 
boundary layer thickness on the asymptotic properties of 
the jet and the ambient medium. 

These models also show, however, that the contact dis¬ 
continuity separating the jet and its surroundings, neces¬ 
sary for maintaining their respective identities, results in 
the likely unphysical vanishing of the density of scatterers 
along that surface of separation. The contact discontinu¬ 
ity also prevents the jet from entraining ambient mate¬ 
rial; since the jet in the two-stream model we considered 
had an inhnite amount of momentum, any entrainment 
or lack thereof is formally inconsequential to the evolu¬ 
tion of the system. However, realistic jets - those with 
finite extent - will almost certainly engulf more mate¬ 
rial as they expand into their surroundings; because the 
total amount of momentum in the system, which is real¬ 
istically finite, must be conserved, that entrainment will 
then result in an overall deceleration of the outflow that 
cannot be captured in the two-stream model. 

In view of these unphysical properties of the two- 
stream treatment - the vanishing of the mass density 
of scatterers along the contact discontinuity and the lack 
of entrainment - we present here an alternative bound¬ 
ary layer scenario to describe the interaction of a rela¬ 
tivistic, radiation-dominated jet with its surroundings. 
This “free-streaming jet” model, which has a well-known 
counterpart in the non-relati vistic, incompressible limit 
(see Chapter 10, Section 12 of iKundu fc Cohenll2008[l . as¬ 
sumes that the jet is injected through a narrow opening 
into a static, homogenous medium, and that far enough 
from that opening the entire system can be modeled as 
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a single fluid that is independent of the details at the 
injection point. By considering the jet and the ambient 
medium as one fluid we obviate the need for a contact 
discontinuity, which we demonstrate allows the density 
of scatterers to remain finite throughout the flow and for 
the jet to entrain material. 

In section 2 of this paper we present the equations of 
radiation hydrodynamics in the viscous limit. Section 
3 uses those equations to analyze the free-streaming jet 
model, and we demonstrate the existence of approximate 
self-similar solutions in the limit that the interaction be¬ 
tween the jet and the ambient medium is concentrated 
in a thin boundary layer. In section 4 we discuss the 
implications of our model and make comparisons to the 
two-stream scenario, and in section 5 we conclude and 
consider the application of the free-streaming jet model 
to super-Eddington TDEs, GRBs, and other astronomi¬ 
cal sources. 


four-velocity I/,y, which giv es (see equation (50) of 
iCoughlin fc Begelmanll20lSl 
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To close the system, we require that the normalization 
of the four-velocity be upheld and that particle flux be 
conserved: 


2. GOVERNING EQUATIONS 

When changes in fluid quantities over the mean free 
path of a photon are small, radiation behaves like an 
effective viscosity and transfers momentum and energy 
between neighboring fluid elements. The precise form 
of the viscosity can be determined by investigating the 
general relativis tic Boltzmann equatio n , which was re¬ 
cently done by ICoughlin &: BegelmanI (|2014bf l for the 
case where Thomson scattering dominates the interac¬ 
tions between the photons and scatterers in the fluid 
rest frame. In this limit, they found that the relativistic 
equations of radiation hydrodynamics for a cold gas (gas 
pressure much less than the gas rest mass density and 
radiation pressure) are (see their equation (49); see also 
Paper I): 


= - 1 , 


p'u^^ 


= 0 . 


( 4 ) 

( 5 ) 


Equations (HD and (|3D ^ (ED constitute six linearly inde¬ 
pendent equations for the six unknowns C/^, e', and p'. 

In addition to the energy density of the radiation, e', 
one can also calculate the number density of photons 
by requiring th at the number flux, be conserved. 
One can show (ICoughlin fc BegelmanI I2014bll that the 
equation = 0 becomes, in the viscous limit, 
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Here the speed of light has been set to one. Creek in¬ 
dices range from 0-3, p' is the fluid rest frame mass 
density of scatterers, e' is the fluid rest frame radiation 
energy density, k is the scattering opacity (in units of cm^ 
g“^), is the metric of the spacetime, is the co¬ 
variant derivative, is the four-velocity of the flow, and 
H'^'^ = U^U’^ +g^’' is the projection tensor. The Einstein 
summation convention has been adopted here, meaning 
that repeated upper and lower indices imply summation. 
This eq uation also shows, in agreement with pr evious 
findings (jBlandford et al.iri985l : lLoeb fc LaoHll992ll . that 
the coefficient of dynamic viscosity, rj, for an optically- 
thick, radiation-dominated gas is 


^ 27 p'k 
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We will also write down the gas energy equa¬ 
tion, obtained by contracting equation (HD with the 


( 6 ) 


where N' is the rest-frame number density of photons. 
Once we solve the equations of radiation hydrodynamics 
0. for the four-velocity of the fluid and the mass density of 
scatterers, we can solve equation ([HD for the number flux 
(.^lof photons. 

The goal of the next two sections is to apply equations 
m and (HD ^ (ED to the boundary layers established be¬ 
tween fast-moving jets and their ambient media. For a 
more thorough discussion of the nature of the equations 
of radiation hydrodynamics in the viscous li mit, we refer 
the reader to (Coughlin fc BegelmanI (|2014bf) . 

3. JETTED BOUNDARY LAYER 

In this section we consider the problem where a nar¬ 
row stream of material is continuously injected into a 
plane-parallel, static, ambient medium, known as the 
free-streaming jet problem. If the Reynolds number of 
the outflow is high, the transition between the stream of 
material and the external environment will be confined 
to a thin layer, permitting the use of a boundary layer 
approximation. 

The basic setup is similar to that of the two-stream 
problem (see Paper I), with the motion of the 2-D, 
plane-parallel injected stream predominantly along the 
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z-direction, the majority of the variation along y, no 
variation or velocity in x, and the point of injection at 
y = z = 0. Now, however, there is no contact discon¬ 
tinuity between the stream and the ambient environ¬ 
ment, meaning that the entire system is considered a 
single fluid. We therefore have less freedom in prescrib¬ 
ing the asymptotic characteristics of the jet and the en¬ 
vironment; however, this configuration permits mixing 
between the two media, allowing the outflow to entrain 
material, which almost certainly occurs in realistic jets. 

We can reduce the complexity of equations (HD and ([3|) 
by assuming that the interaction between the jet and 
the ambient medium takes place over a thin boundary 
layer of thickness Sy, thin in the sense that S = SyjSz 
is a small number when 6z is a typical length along the 
jet. In this case, S scales in an identical fashion to that 
derived for the two-stream boundary layer, which can 
be determined by comparing leading-order terms in the 
boundary layer thickness to the inviscid terms in the gas 
energy equation. Setting the left-hand side of equation 
(ini) to ^ e'TjVj/z, we find, as in Paper I, 


S' = In -I- \/l + = arcsinh(ru 2 ). (14) 

The first of these shows that the radiation energy den¬ 
sity, and consequently the pressure, is constant across 
the boundary layer. In the inviscid limit, equation (I12p 
can be transformed to give 


c/^v^ 


r 



= 0 , 


(15) 


which, as we mentioned, is the relativistic generalization 
of the Bernoulli equation. We will also assume that the 
ambient energy density is independent of z, i.e., that 
e'(z) = e'g, with e'g a constant. 

Because it will be convenient, we will change variables 
from y to T, where r is given by 


rv 

p\y,z)Kdy, (16) 

^0 


p' KZ Pj Vj 

Here Vj and P^- = (1 — are the jet velocity and 

Lorentz factor, respectively, measured at some character¬ 
istic length along the jet axis zg, and p' is a characteristic 
density of scatterers throughout the outflow. Due to the 
fact that the boundary layer thickness is the same, the 
boundary layer equations governing the outflow are iden¬ 
tical to those found in paper I, which can be compactly 
written as 






3 ^ 27 i9y \_p'k dy J ’ 

(8) 

p'U>^ = 0. (9) 


which is related to the optical depth across the boundary 
layer as measured from the axis. In terms of this variable, 
equations m and (fT^ become, respectively. 
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where we have used the assumption that e'{z) = e'g. It 
should also be noted that the y-component of the con¬ 
vective derivative is now with respect to r, not y, i.e.. 


= (19) 

In the ensuing section we will seek self-similar solutions 
to these equations. 


In paper I we dealt with the v = y and v = z compo¬ 
nents of equation ([8]) and its contraction with Uv - the 
gas energy equation ©■ For the free-streaming jet prob¬ 
lem, however, it will be more convenient to deal with 
the u = y component of equation ([5]), the gas energy 
equation, and the contraction of equation (|5]) with the 
projection tensor which, as we will see, is the rel¬ 
ativistic, viscous counterpart of the Bernoulli equation. 
Performing a few manipulations, we find that these equa¬ 
tions become, respectively. 


dy 


= 0 , 


e'U'' 
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3.1. Self-similar approximation 

We can immediately solve the continuity equation m 
by introducing the stream function if via 


PUz = 


d-tp 


p'nTvy = 


dtp 


( 20 ) 

( 21 ) 


where the factor of k ensures that tp remains dimension¬ 
less. In paper I we showed that there exist self-similar 
solutions for tp, the velocity, the comoving density of 
scatterers and comoving density of photons in terms of 
the variable a = y/Sy, where 6y ^ ^/z. One difference 
between the two-stream problem and this type of out¬ 
flow, however, is that we expect the jet to expand into 
the ambient medium, entraining material in the process. 
Therefore, as we look farther along the z-direction, the 
amount of inertia - predominantly in the form of radia¬ 
tion for the systems that we are considering - contained 
in the flow will increase. Owing to the conservation of 
momentum, the z-component of the velocity should thus 
be a decreasing function of z. 


where 
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In light of this observation, we will assume that the 
z-coniponent of the four-velocity scales as 


= TjVj 





( 22 ) 


where zq is a characteristic length scale in the z-direction, 
TO is a positive constant, for brevity we dehned = 
TjVj{z/zo)~'^, and / is a function of the self-similar vari¬ 
able f ~ tI8t\ 5t is the characteristic boundary layer 
thickness in terms of r, which, from equations © and 
dm) in the previous subsection, is given by 5t ~ p'kSi/, 
or 


tained in the outflow be conserved. Integrating 
from y = —oo to oo, using equations (1^^ and ((Ml) , tak¬ 
ing the ultrarelativistic limit and changing variables from 
y to we hnd 



Since the integral in this equation is a constant that is 
greater than zero, we hnd that to = 1/4 if the momentum 
hux is conserved. 

With this value of to, the z-component of the four- 
velocity scales as 




p KZ 


(23) 


We are most interested in the behavior of the prop¬ 
erties of the outhow when the velocities are relativistic, 
as this is the limit that is most applicable to sources of 
astronomical interest, and in order for the self-similar na¬ 
ture of our solutions to be upheld, the dependence of our 
self-similar functions on the bulk properties of the out¬ 
how, such as the z-dependence of the jet Lorentz factor 
along the axis, should be minimal. A fully self-similar so¬ 
lution is likely impossible here, as the speed of light plays 
a role in setting a hnite scale factor for our solutions and 
becomes problematic when we try to connect the ultra- 
and non-relativistic regions of the outhow. We will show, 
however, that the non-self-similarity of our solutions only 
affects a small region of the outhow. 

In paper I, we showed that the self-similarity of the co¬ 
moving density of scatterers, p', was approximately sat- 
ished, i.e., we could hnd solutions with p' ^ g{0- For the 
case at hand, then, this assumption implies that g is in¬ 
dependent of C. This means, however, that the observer- 
frame density of scatterers, given by p = Fp' ~ Cp, can 
be made arbitrarily large (as g is independent of in the 
self-similar limit), which is a nonsensical result. 

Motivated by this reasoning, we conclude that the co¬ 
moving density of scatterers is unlikely to be independent 
of C- On the contrary, a more reasonable approximation 
is that the observer-freime density of scatterers varies 
self-similarly. Investigating equations dni) and m , we 
see that the p'-dependent quantity that enters both is 
the combination 1 -|-4e'/(3p'). Our expectation that the 
observer-frame density varies self-similarly then prompts 
the assumption 


l + ~ = ^gTg{0, (24) 

3 p' 3 

where p = e' (Pq, Pq being the density of scatterers in the 
ambient medium. (Note that we could have simply let 
p'r = g{^), but equation (|M)) . which is merely a change 
of variables from the initial assignment p'F = p(^), will 
allow equations dni) and m to be written in a more 
compact form.) 

We can determine the value of to, which controls how 
rapidly the flow decelerates due to entrainment, by inte¬ 
grating the momentum flux, P ~ e'F^v^ for a radiation- 
dominated system, over the entire boundary layer, and 
requiring that the result be independent of z. This is 
equivalent to requiring that the total momentum con- 


rvz = 




and the self-similar variable ^ is given by 


(26) 


^ = T- 


(- 

2pgKZo V-^O 




-3/4 


(27) 


We then hnd from equation (|M)) that the stream function 
must satisfy 


^ = (28) 

Equation (IMI) then gives the p-component of the four- 
velocity: 


p'nTvy = — 





. (29) 

where, both here and in future equations, a subscripted 
^ on the functions / and g denotes differentiation with 
respect to ^ and the number of subscripts indicates the 
number of derivatives, i.e., /j = df/d^, = d'^f/d^^, 
etc. 

Inserting equations ((Ml) . (IMl) and (l29l) into equations 
(HZD and m, changing variables from r to ^ and per¬ 
forming a bit of algebra, we hnd that they become, re¬ 
spectively, 




- - 


d^s 

ae ■ 


(30) 

(31) 


Now, the solutions we seek for / and g should de¬ 
pend only on as otherwise our assumption of the self¬ 
similarity of those functions breaks down. Because of the 
complicated dependence of S' on Fuz, we see that the so¬ 
lutions will not be self-similar for arbitrary Fj. However, 
in the ultrarelativistic limit, for which F ~ Fu^, we see 
that S ~ In(Fu^), and it is apparent that the solutions 
are indeed self-similar, i.e., the C dependence of equations 
((5D1) and ((Ml) drops out. 

As these equations are third order in / and hrst order 
in g, we require four boundary conditions to solve them 
numerically. Since Fuz is given by equation (1261) and it is 
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assumed that Vvz = TjVj{z /along the axis, our 
first boundary condition is given by /j(0) = 1. For the 
second, we note that the jet axis is the streamline along 
which the y-component of the velocity is zero. From 
equation (1^^ . then, we find /(O) =0 (we will see that 
d^jdz = 0 at ^ = 0). The third boundary condition is 
obtained by noting that the 2 -component of the velocity 
should approach zero as we proceed into the ambient 
medium, so that, from equation (l26l) . we find /^(oo) = 0 . 
Finally, the density of scatterers should approach that of 
the ambient medium far from the jet center. Equation 
dm) then gives g{oo) = f -|- 3/(4^). 

As we noted, the truly self-similar limit of equations 
((5D1) and (l?T|) is obtained by setting S = In(Fuz), but we 
encounter an issue with this scaling when we consider the 
boundary conditions on our flow in the ambient medium. 
In particular, we expect the outflow velocity to approach 
zero for ^ ^ 1 ; when the velocity becomes subrelativistic, 
however, ln(F?;z) —>■ —oo, when in actuality the Vz 
I limit of equation (THl) is 5" ~ ln(l -|- Fuz) ~ Tvz — 
0. Thus, we see that taking the ultrarelativistic limit of 
equations dSOl) and m will not result in the solutions 
matching the correct boundary conditions far from the 
jet center. 

To correct this problem, we will allow the function S to 
take on its full form and treat C as a constant, which will 
break the self-similarity of our solutions. However, since 
our boundary conditions do not depend on the value of 
C, the functions / and g themselves should be largely in¬ 
dependent of that parameter. We therefore expect the 
assumption of self-similarity to be upheld in the ultrarel¬ 
ativistic <C 1 ) and the non-relativistic ^ 1 ) limits of 
our solutions, with small deviations from self-similarity 
in the trans-relativistic regime of the outflow. 



Fig. 1. — The behavior of /^, which is the normalized 2- 
component of the four-velocity, with ^ for C = 4 (blue, solid curve), 
(^ = 10 (purple, dot-dashed curve), and <^ = 25 (red, dashed curve). 
As expected, the curves are all coincident when ^ 1 and ^ 1, 

with a non-self-similar transition (one that depends on (^) in be¬ 
tween those two limits. 



Fig. 2. — The function g, which is approximately the inverse of 
the lab-frame density, plotted with respect to ^ for the same set 
of ^ chosen in Figure ^ As was true for /^, g is approximately 
self-similar close to and far from the jet, with the deviation from 
self-similarity in the trans-relativistic region being apparent but 
small. 

relativistic, the deviations from self-similarity are appar¬ 
ent, though they remain small. Keeping the full form of 
S in equations ([3Q|l and ([3T|) therefore preserves well the 
self-similarity of our solutions and provides a reasonable 
interpolation between the relativistic and non-relativistic 
regions of the flow. 

The solution for g approaches infinity as we near the 
axis, which can be seen by investigating the small-^ be¬ 
havior of equation m- Specifically, letting / ~ ^, 
~ 1, F ~ C, and Uz — 1, equation (1311) can be written 
as 


^5? + 2^ — (32) 

where is the second derivative of / evaluated at 
^ = 0. This equation can be integrated, and we find 

9--^- (33) 

C being a constant of integration. Since it is roughly 
proportional to l/p^ the value of g cannot be negative, 
meaning that C must be greater than zero. Therefore, 
the asymptotic, ^ <C 1 behavior of g is 5 ^ I/a/J- 
In addition to the mass density of scatterers and the 
velocity, we can also calculate the number density of pho¬ 
tons, N', throughout the boundary layer. As we showed 
in Paper I, the equation of photon number conservation 
([S]) becomes, to lowest order in the boundary layer thick¬ 
ness. 


N'U^ 


19/1 dN'\ 
idy\p'K dy )' 


(34) 


Performing a few manipulations, this equation becomes 


To exemplify this point. Figures [T] and [5] show, respec¬ 
tively, the variation of - the normalized 2 -component 
of the four-velocity - and g - approximately the inverse 
of the observer-frame density - in terms of the parameter 
^ for C = 4,10, and 25. As we anticipated, the functions 
for different / are indistinguishable for ^ <C 1 and ^ ^ 1 , 
meaning that the self-similarity is nearly exact in those 
regions. In between those limits, where the flow is trans- 





nd^N' 
3 9r2 ■ 


(35) 


As was true for the density of scatterers, we expect the 
observer frame number density of photons to vary ap¬ 
proximately self-similarly, which will be true if we let 


N' 

P' 




( 36 ) 
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where Nq is the number density of photons in the ambient 
medium and h is a dimensionless function. Inserting this 
ansatz into equation (l3^ gives 





(37) 


For the boundary conditions on h, we first require that 
the number density of photons approach that of the am¬ 
bient medium in the ^ ^ 1 limit. Equation (I36p then 
gives h{oo) = 1. For the second condition, return to 
equation (I34|l . integrate both sides from y = —oo to oo, 
and require that the derivative of N vanish in both of 
those limits. Doing so yields 


dz 



N'Tv^dy 


-N'TVy 


OO 


— OO 


(38) 


The right-hand side can be determined by returning to 
the continuity equation, integrating from y = —oo to 
00 and performing a few manipulations to show that 
r'*^oo|_^ = —foe, where /oo is the function / evaluated 
at infinity. Using the definition of N' in terms of h, we 
find that equation (IMll becomes 


h f(_di = 2 f^ 


(39) 


which serves as our second boundary condition on h. 
This integral states that the increase in the number flux 
of photons occurs at a rate provided by the influx of ma¬ 
terial at infinity. 

The equations derived in this section were all written 
in terms of the variable which is itself a function of 
g via equations (HID and (HZD- To write the solutions 
in terms of the physical parameter y, we can return to 
equation (I16|l . differentiate both sides with respect to y, 
rearrange the resulting equation and integrate to yield 




di = y 




-3/4 


= a. 


(40) 


where ^ is a dummy variable of integration. Once we 
calculate the functions g{^) and this relation can 

be integrated and solved numerically to yield ^(a). This 
expression also shows that 


rg 


Afijdz 4 z\J, v '1 + C2(/5)2 

which we can use in equation (1291) to give 


d^ — 3a I, (41) 


p'nTvy = 

4 V 9 z \Tg-^ 



e{k?9 

x/TTcW 




Equation (ITip also confirms that d^/dz = 0 when / = 
0 , which we used in order to determine the boundary 
condition /(O) = 0. 


3.2. Solutions 

In this section we plot solutions for the outflow velocity, 
the density of scatterers and the density of photons for 
various values of g and /. As was mentioned in the previ¬ 
ous subsection, the physical self-similar variable against 
which we would like to plot our solutions is given by a. 
However, as is apparent from equation (14011 . the defini¬ 
tion of a depends on /. Therefore, if we are comparing, 
for example, the outflow velocity of two systems with dif¬ 
fering we must incorporate the / dependence in a so 
that the range of physical space that we consider for each 
solution is the same. For this reason, in this section we 
will plot our solutions as functions of the variable 



^ p'qK f z 

2 zo V-^o 


- 1/2 


(43) 



a 

Fig. 3. — The normalized ^-component of the four-velocity (/^) 
for /r = 1 and (^ = 4, 10, and 25 (the solid, blue curve, the dot- 
dashed, purple curve, and the dashed, red curve, respectively), 
which, for z ~ 20 , correspond to = 4,10, and 25. We see that 
the width of the boundary layer is nearly unchanged as we alter 
the value of 



a 

Fig. 4.— The normalized fluid-frame density of scatterers for 
the same set of parameters used in Figure 0 For all solutions the 
number density of scatterers approaches zero as we near the center 
of the jet. We see that the average number density of scatterers 
within the boundary layer is lower for larger Lorentz factors. 

Figure [3] shows the solution for the normalized z- 
component of the four-velocity, for / = 4, 10, and 25. 
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a 

Fig. 5. — The normalized number density of photons for the 
same set of C, used in Figure [3| The photon number density closely 
follows that of the scatterers. 



a 

Fig. 6. — The ^-component of the three-velocity normalized by 
-\/2/(9 pqK2:o) (^e equation (142 l l) for the same set of parameters 
used in Figure [3] For positive a, each solution is initially positive 
and reaches a relative maximum before approaching a negative 
constant, which shows that the flow expands outwards near the 
center of the jet and entrains ambient material far from the axis. 

Since C = ^jVjiz/, these values of ( scale approx¬ 
imately linearly with the Lorentz factor until z ^ zq. 
The outflow velocity is maximized at the origin and de¬ 
cays as we move farther into the ambient medium. We 
see that the width of the boundary layer, loosely dehned 
as the value of a at which /j is some fraction of its central 
value, is nearly unchanged as we modify The average 
value of the normalized z-component of the velocity is 
also slightly larger for larger 

Figure 13] demonstrates how the normalized, fluid-frame 
mass density of scatterers varies as we traverse the 
boundary layer for the same set of parameters chosen 
in Figure |3l Since g approaches infinity as we near the 
origin, the mass density of scatterers, related to g by 
equation (IMl) . equals zero at the origin for all of the so¬ 
lutions, meaning that the center of the jet is evacuated of 
massive particles. The average comoving density of scat¬ 
terers across the boundary layer is also lower for larger C, 
in accordance with equation ((24ll . In Figure [5] we plot the 
normalized number density of photons for the same set of 
(■. It is evident that the density of photons closely follows 
the density of scatterers throughout the boundary layer. 
We see, however, that the density of photons stays above 
and below the density of scatterers as we move toward 


and away from the center of the jet, respectively. The 
photon density also remains finite at the center of the 
jet. 

The y-component of the three-velocity normalized by 
yj2/{9pQKZo) ~ (5 is illustrated in Figure [HI for the same 
set of C used in Figure |31 For ^ > 0, each solution ini¬ 
tially has a positive Uy, which shows that the jet mate¬ 
rial expands away from the axis. The y-component of 
the three-velocity then reaches a relative maximum, one 
which increases slightly for larger F^ , before approach¬ 
ing a negative, constant value. This behavior is then 
inverted for negative Because the transverse velocity 
approaches a negative constant for ^ > 0 and a positive 
constant for ^ < 0, we see that the jet entrains material 
from the ambient medium. 



a 

Fig. 7. — The function /^, which is the normalized .^-component 
of the four-velocity, for ^ = 10 and n = e' j= 0.1, 1, and 10, 
which correspond to the blue, solid curve, the purple, dot-dashed 
curve, and the red, dashed curve, respectively. Increasing the value 
of /i., we see, has little effect on the solution, while decreasing p 
drastically widens the boundary layer. 



a 

Fig. 8. — The normalized density of scatterers for the same val¬ 
ues of p chosen in Figure^! The mean value of the density decreases 
as p decreases. 

Figures 13 -(TOj illustrate how our solutions depend on 
fi. As is apparent, changing the value of p, does not dras¬ 
tically alter the qualitative aspects of the functions. We 
do see, however, that decreasing p from 1 to 0.1 results 
in a large increase in the boundary layer thickness; on 
the contrary, changing p from 1 to 10 results in only a 
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(44) 



a 

Fig. 9. — The normalized number density of photons for the 
same set of /i chosen in Figure [7| This figure demonstrates, as 
we saw in Figure O that the density of photons tracks that of the 
scatterers. 



a 

Fig. 10. — The y-component of the three-velocity, normalized 
by ^y2/{9pQKZQ) ~ 5, for the same set of p chosen in Figure [71 We 
see that the relative maximum increases for smaller p. 

slight narrowing of its thickness. It is also evident that a 
smaller compared to 1 results in a lower average value 
of the density throughout the boundary layer and a larger 
peak in the transverse velocity Vy. 

4. DISCUSSION 

Figures[3]-[n]demonstrate that the width of the bound¬ 
ary layer is nearly independent of which is the re¬ 
sult of the scaling of our boundary layer thickness (5. 
Specifically, note from equation m that S is given by 
6 ^ l!^p'Tj. Since our ansatz posited that p'T ~ g, 
with g a dimensionless function of order unity, and 
r — Tj ~ C for ^ the boundary layer thickness 

is roughly independent of C 

Figures [3-[To] illustrate that a value of p = 0.1 causes 
the boundary layer thickness to increase dramatically 
compared to p = 1 , while setting p — 10 causes only 
a slight narrowing of the width compared to p = 1. This 
dependence is due to the fact that the viscous heating, 
which increases the specific entropy by decreasing the 
density of scatterers (at fixed pressure), is most efficient 
when the flow is compressible, as is evident from the 
gas energy equation m- The fluid only becomes com¬ 
pressible, however, when the flow is sup ersonic, and we 
can show (jCoughlin fc Begelmanll^014bl) that the sound 
speed of a radiation-dominated gas is 


l + 4p/3 

When p <C 1, the sound speed reduces to Cs ^ y/p, 
and the location at which the flow becomes transsonic 
extends farther into the ambient medium, widening the 
boundary layer. Conversely, when p ^ 1, the sound 
speed approaches a constant ~ 1/v^, which results in 
only a slight narrowing of the boundary layer. 

Equation (|10p shows that the energy density of the ra¬ 
diation, and hence the pressure, is constant across the 
jet, which is a statement of the causal connectedness of 
the boundary layer. In order for this equation to remain 
valid, then, we require that the transverse sound crossing 
time over the boundary layer thickness 6y be less than the 
time it takes the fluid to traverse the distance Sz. Since 
the transverse sound speed is given c_l = Cg/T, where Cg 
is given by equation (|44l) , this requirement yields the in¬ 
equality S ^ Cs/iTj Vi). Once this inequality is no longer 
satisfied, equation (| 10 p does not hold, and we must in¬ 
clude more terms in equations (uni) - dm) to account for 
the gradients in the energy density of the radiation. 

As we noted in Section 3, the density approaches zero 
as we near the axis of the jet. Physically, this effect 
arises from the fact that, for z zq, the Lorentz factor 
grows unbounded and the boundary layer thickness goes 
to zero. Therefore, the center of the jet originates along a 
curve of infinite shear and, consequently, infinite entropy. 
Since the specific entropy scales as s' oc e'/(p')^/^ and e' 
is constant across the boundary layer, we see that the 
density of scatterers must equal zero at the center of the 
jet. 

From Figures |5| and (TUI we see that the j/-component 
of the velocity initially causes the outflow to expand into 
its surroundings. When |^| becomes large, however, the 
directionality of Vy reverses toward the jet. The outflow 
therefore entrains material from the ambient medium not 
only by expanding in the transverse direction, but also 
by dragging material in from the environment. 

Interestingly, the y-component of the velocity ap¬ 
proaches a non-zero, constant value as we move into the 
external medium. This behavior arises from the fact that 
the jet is removing material from the system at a rate 
M ~ p'Tvzdy ~ showing that amount of mass ex¬ 
cavated from the envelope increases as we move along 
z. Therefore, in order to maintain a steady-state, we re¬ 
quire a constant influx of material at infinity that can 
resupply the amount lost due to the jet. This inter¬ 
pretation is substantiated by integrating the continuity 
equation Cl from y = —oo to y = oo, which shows 
that the y-component of the velocity at infinity scales as 
Vy,oo ^ —/oo/2. The factor of 1/2 arises from the fact 
that the mass loss rate scales as M ^ p'Tvzdy oc 
therefore, if one could decrease the mass loss rate to one 
that was constant in z, then the y-component of the ve¬ 
locity would vanish at infinity. Likewise, if one could 
create a scenario in which M scaled as a negative power 
of z, then the value of Vy would maintain an efflux of 
material at infinity to keep the system from amassing 
inertia towards the center of the jet (thus violating the 
steady-state assumption). 

The free-streaming jet solutions analyzed in this paper 
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are only valid when the Lorentz factor is relativistic, i.e., 
when C > 1- When C<Cl,r~l, p'~/5 q, and we can 
show that conservation of momentum along the z-axis 
results in the scaling Vj ^ S ^ which agrees 

with the results of the incompressible , non-relativistic 
theor y (in Chapter 10, Section 12 of iKundu fc CohenI 
I2008L see their discussion at the top of page 383). There¬ 
fore, once z zq, the outflow will undergo a non-self¬ 
similar transition from the solution presented here to its 
non-relativistic counterpart. 

The free-streaming jet model has observational conse¬ 
quences. For example, the decrease in the density along 
the jet axis means that observers looking down the barrel 
of the jet see farther into the outflow. Because T oc z~^^^ 
and e = T^e', where e is the lab-frame radiation energy 
density, this means that those observers see a higher en¬ 
ergy density of photons. Also, those same observers, us¬ 
ing the fact that the lab-frame radiation number density 
is N ^ TN', see an energy per photon oi e/N ^ Te'/N'. 
Therefore, not only do the on-axis observers see a more 
Lorentz-boosted spectrum because they can see deeper 
into the outflow, their spectrum is also hardened from 
the fact that N' is minimized along the axis of the jet 
(where j/ = 0 in the figures of Section 3). 

4.1. Comparative notes between the two-stream model 
(Paper I) and the free-streaming jet model 

Paper 1 analyzed how the presence of a radiation 
pressure-supported envelope affected the propagation of 
a radiation-dominated, relativistic jet through the two- 
stream approximation. This approximation treats the 
jet and its surroundings as semi-infinite, separate fluids, 
a contact discontinuity serving as the surface of separa¬ 
tion between the two. In contrast, the free-streaming jet 
solution presented here treats the entire system - the jet 
and its surroundings - as a single fluid. There is thus no 
formal distinction between “jet material” and “ambient 
material,” meaning that no contact discontinuity exists 
in the system. 

By comparing Figure 1 of Paper I and Figure |3] of the 
previous section, we see that the z-component of the four- 
velocity at a fixed z behaves similarly between the two 
models. Namely, the solution starts at some “jet” veloc¬ 
ity, which corresponds to y = oo in the two-stream model 
and to y = 0 in the jet model, and smoothly transitions 
to a velocity of approximately zero over the extent of a 
few boundary layer thicknesses. However, the full spatial 
dependence of the velocity, one that includes variation in 
the z-direction, differs drastically between the two mod¬ 
els: while the two-stream jet maintains a constant Fj 
along z, we found here that the z-component of the four- 
velocity scales as rz;^ oc implying that the over¬ 

all velocity of the jet slows as we look farther down z. 
This behavior arises from the fact that the free-streaming 
jet can entrain material, this entrainment causing an in¬ 
crease in the inertia contained in the outflow and a re¬ 
sultant decrease in its velocity. 

The general behavior of the number densities of scat- 
terers and photons between the models is also similar, 
which can be understood by comparing Figures 2 and 3 
of Paper I to Figures |T] and [S] of Section 3, respectively. 
In particular, both densities decrease within the bound¬ 
ary layer separating the outflow and its surroundings. 


asymptotically approaching their jet values as y —> cx) in 
the two-stream model and as y —>■ 0 in the free-streaming 
jet model. Likewise, each approaches its ambient value 
as y —>■ —oo in the two-stream model and as y —>■ ±cx) in 
the free-streaming jet model. Although the gross prop¬ 
erties of both are similar, one striking difference arises, 
however, in the behavior of the scatterer density: in the 
two-stream model, the existence of the contact disconti¬ 
nuity causes the function g, and hence the density p', to 
vanish within the boundary layer (it vanishes specifically 
at y = 0 - the location of the contact discontinuity). 
Conversely, p' remains non-zero throughout the bound¬ 
ary layer that connects the jet and its surroundings in 
the model presented here, and only as we near the cen¬ 
ter of the jet does the density of scatterers go to zero. 
The presence of a contact discontinuity thus results in 
the likely non-physical vanishing of the massive particles 
within the boundary layer. 

Finally, the scaling of the boundary layer thickness it¬ 
self differs between the two models. In the two-stream 
case, we found that 5 oc z^/^ (see equation (15) of paper 
I). Therefore, as one looks farther down the z-axis, the 
boundary layer that develops between the jet and the 
ambient environment extends into both media at a rate 
proportional to z^^^. Contrarily, equation (SOI) shows 
that the free-streaming jet boundary layer expands into 
its surroundings as oc z^/'^ (though this transitions to z^/^ 
in the non-relativistic limit; see the dis cussion at the top 
of pa ge 383 in Chapter 10, Section 12 of iKundu fc CohenI 
I2008f) . The free-streaming jet boundary layer thus ex¬ 
pands more rapidly than does the two-stream solution. 

5. SUMMARY AND CONCLUSIONS 

Employing the equations of radiation hydrodynamics 
in the viscous limit, which are applicable as long as 
changes in fluid quantities are small over the mean free 
path of a photon, we analyzed the dynamics of a relativis¬ 
tic, free-streaming jet under the boundary layer approx¬ 
imation. This approximation, which should be upheld in 
jets with transverse optical depths substantially greater 
than one, states that variations in the properties of the 
outflow are confined to a thin layer of width 6, and it 
allowed us to transform the full set of equations o-® 
into a set of greatly simplified boundary layer equations 

(0 - ini)- 

Perhaps the biggest difference between the two-stream 
solutions, presented and analyzed in paper I, and the 
free-streaming jet solutions presented here is in the dis¬ 
tinction between the jet and the ambient medium. In 
the former, the two are considered as distinct, interact¬ 
ing entities, which allows one to specify separately their 
asymptotic properties. The latter approach, on the other 
hand, considers the whole configuration as a single fluid. 

Because one has more freedom in specifying the prop¬ 
erties of the outflow, the two-stream solution has the 
added benefit of being able to treat scenarios in which 
the properties of the jet and the ambient medium dif¬ 
fer significantly. However, maintaining the distinction 
between the two media necessitates the existence of a 
surface of contact between the two across which no fluid 
can flow, meaning that the jet cannot entrain ambient 
material. Furthermore, as was demonstrated in Paper I, 
this boundary condition results in the density formally 
vanishing at the interface, which is likely non-physical. 
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By treating the outflow and the environment as one 
fluid, we demonstrated that entrainment does occur in 
the radiation-viscous, free-streaming jet solution, which 
causes the z-component of the four-velocity of the jet to 


slow as ^ (this power-law, however, may differ if 

the adopted symmetry is azimuthal as opposed to planar; 
see below). These solutions also show that the comoving 
densities of scatterers and photons decrease dramatically 
within the boundary layer, with the density of scatter¬ 
ers approaching zero as one nears the center of the jet. 
Therefore, because observers that look “down the bar¬ 
rel of the jet” can see farther into the outflow, they see 
a more Lorentz-boosted energy density than those that 
view the outflow off-axis. They also observe a higher 
energy per photon, given by e/N ~ e'To/N', both be¬ 
cause N' is minimized along the axis and because To, 
the observed Lorentz factor, is larger. Such features are 
in qualitative agreement with the event Swift J1644+57, 
where such an obser ver orientation is inv oked to explain 
the X-ray emission (jZauderer et al.l[2011ll . 

In addition to super-Eddington TDEs and long GRBs 
(as we ll as the relativ ely new class of “ultra-long” 
GRBs; ILevan et al.ll201'3l - the two applications consid¬ 
ered in the Introduction - the free-streaming jet model 
could be applied to other sources. As mentioned in 
Paper I, this model may also be relevant to micro¬ 
quasars ([Fender et al.l 120 0411 fparticular ly those that fall 
in the class o f ULX’s ; King et_aD [2001D such as the ob- 
iect SS 433 (iFabrikal 120041 : iBegelman et all l2006al: see 
lArav fc BegelmanI (|1993f) for an application of the two- 
stream model to this source). Additionally, a jetted 
quasi-star - a protogalactic gas c loud supported by a su- 
percritically accreting bla ck hole ([Begelman et al.ll2006bL 
120081 iCzernv et al.ll201^ - provides another situation in 
which a radiation-dominated jet propagates alongside a 
radiation pressure-supported envelope. 

The free-streaming jet solution presented here is lim¬ 
ited to describing plane-parallel, two-dimensional sys¬ 


tems. One consequence of this assumption is that the 
entrainment of material, which slows down the jet as 
r oc happens effectively in one dimension. If, on 

the other hand, one imposed azimuthal symmetry and 
described the system in terms of spherical coordinates 
(r, 0, (()), which is likely more relevant for realistic jetted 
systems, the entrainment would occur in two dimensions. 
Other conditions being equal, the slowing of the jet along 
the axis would then be more pronounced, i.e., the Lorentz 
factor would scale as T oc r~P with p > 1/4. 

The other main assumption of the free-streaming jet 
model presented here is that the energy density of the 
ambient medium is independent of z, i.e., e'(z) = eg. If 
a pressure gradient were present, this force would tend to 
accelerate (or, in principle, decelerate) the jet material, 
offsetting the power-law scaling T cx In fact, one 

can imagine that if the pressure gradient were strong 
enough, it could indeed reverse the overall slowing of the 
jet and cause the outflow to accelerate. If the energy 
density scaled as e' = eg(z/zo)“'^, which is relevant for 
super-Edd ington TDEs as long as the con stant q satisfies 
q > 3/2 (I Coughlin &: BegelmanI [Ml4all . this situation 
could be actualized near the launch point of the jet. 

In an ensuing paper, we plan to compare more quan¬ 
titatively the predictions made by the free-streaming 
jet model presented here and the observations of Swift 
J1644-I-57. We will also extend our analysis to incor¬ 
porate a spherical geometry - as azimuthal symmetry 
is almost certainly more relevant for this and other sys¬ 
tems than the planar symmetry adopted here - as well as 
radially-dependent ambient pressure and density profiles. 


This work was supported in part by NASA Astro¬ 
physics Theory Program grant NNX14AB37G, NSF 
grant AST-1411879, and NASA’s Fermi Guest Investi¬ 
gator Program. 
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